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Abstract 

Given measurement data, a nominal model and a linear fractional transformation uncer- 
tainty structure with an allowance on unknown but bounded exogenous disturbances, easily 
computable tests for the existence of a model validating uncertainty set are given. Under 
mild conditions, these tests are necessary and sufficient for the case of complex, nonre- 
peated, block-diagonal structure. For the more general case which includes repeated and/or 
real scalar uncertainties, the tests are only necessary but become sufficient if a collinearity 
condition is also satisfied. With the satisfaction of these tests, it is shown that a parame- 
terization of all model validating sets of plant models is possible. The new parameterization 
is used as a basis for a systematic way to construct or perform uncertainty tradeoff with 
model validating uncertainty sets which have specific linear fractional transformation struc- 
ture for use in robust control design and analysis. An illustrative example which includes a 
comparison of candidate model validating sets is given. 

1 Introduction 

In applying multivariable robust control analysis and synthesis techniques to linear, time- 
invariant systems, a particular set of plant models described by a nominal model, uncertainty 
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structure, and norm bounds on the model uncertainty and exogenous inputs is required a 
priori [1, 2], In a typical setting, a controls engineer selects or develops a best possible model 
either from first principles and/or from system identification or parameter estimation. In 
situations where the physical conditions are not accurately known nor reliable, the engineer 
selects a model uncertainty structure around a best possible nominal model in order to rep- 
resent a set of plant models so that robustness of the system can be analyzed and optimized 
through feedback control. In many instances where measurements are available, the engineer 
can check to see if the given set of plants is consistent with the available measurement data. 
Specifically [3] -[4], given a robust control design model, does a plant exist within this set 
which will reproduce a given input and output measurements, i.e. , is the given set of plants 
model validating (or more accurately cannot be invalidated) with respect to available data? 

There exists a wealth of literature related to model validation, uncertainty modeling and 
identification for control (see, for example [5]). Although there have been several successful 
methodology validation tests reported recently as in [6]- [10] , it appears that there is not a 
systematic methodology that can be demonstrated to work for a significant class of problems. 
Indeed, this study addresses a variant model validation problem originally formulated in 
[3]. In [3], a smallest level of exogenous disturbance that satisfies model validation given a 
fixed level of model uncertainty bound is sought whereas in this paper, we seek a smallest 
model uncertainty bound that satisfies model validation given an allowance on the exogenous 
disturbance. Furthermore, exogenous noise signals of the type Bi ^ is assumed in [3] while 
we assume that the exogenous noise is a colored white noise signal whose power spectrum is 
a known constant over a wideband frequency, i.e. belongs to B^°°\ 

The model validation question has a binary outcome which can be viewed as a check 
on two properties. First, the postulated parameteric and nonparametric uncertainties and 
their corresponding interconnections leading to a Linear Fractional Transformation (LFT) 
structure is required to be sufficiently rich to admit a perturbed model that will faithfully 
reproduce the measured input/output data. In the seminal work in [3], this property is 
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referred to as a feasible set. Second, assuming feasibility, the postulated norm bound on the 
uncertainties must be sufficiently large to admit a model validating plant. The viewpoint 
taken in this paper with regards to model validation is to focus on the first property, namely 
whether an a priori given LFT structure for an uncertainty model can lead to a model 
validating set given some allowance in exogenous disturbance. The rationales for the above 
viewpoint are: (i) feasibility question is simpler to address, (ii) if feasibility is not satisfied, 
there is no point in checking the second property, (iii) once feasibility is satisfied, a model 
validating set can almost always be constructed, using the results in this paper. The above 
viewpoint suggests a slight adjustment in the model validation paradigm, namely, it is not 
of prime significance whether a particular set satisfies model validation since if feasibility is 
satisfied, model validating sets are highly non-unique. 

For brevity, we highlight only closely related earlier work to our approach. In an early 
attempt to obtain a simpler problem formulation and solution than the approach taken by [3] 
to [6], [8] considered a special case where all exogenous inputs are either known or are very 
small and occur only at the output. This attempt was subsequently applied to an experi- 
mental testbed with encouraging results [9]. Although the approach taken in [8], [9] appears 
to work reasonably well for problems with an arbitrary number of structured full complex 
blocks only, it became clear through applications that problems with parametric (and often 
real repeated) uncertainties gave unsatisfactory results. This was somewhat expected be- 
cause the additional structure in the repeated scalar uncertainties was not incorporated in 
the original problem formulation. Hence, the most recent work reported in [11] extends the 
previous approach to include repeated and/or real scalar parametric uncertainties along with 
an arbitrary number of full complex blocks. This paper provides the detailed proofs omitted 
in [11], expands the uncertainty set to include unknown but bounded output noise and/or 
disturbance to the plant, and ultimately parameterize all model validating sets for a fixed 
LFT structure. We also discuss and provide a means to cope with the issue of determin- 
ing non-parametric uncertainty bounds in the presence of constant parametric uncertainties. 
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The formulation in this paper is for a closed loop system and the open loop configuration is 
obtained as a special case. 

In Section 2, a problem definition is given whereby uncertainty bounds are viewed as 
bounds on fictitious uncertainty signals which satisfy P — A transmission conditions while 
resulting in zero output errors. In Section 3 we derive existence conditions for a model vali- 
dating set followed by a parameterization of all model validating sets and associated signals. 
In Section 4, we outline possibilities in utilizing the parameterizations given in Section 3. A 
smallest set formulation is discussed in which allowances are given to parametric uncertain- 
ties similar to the treatment of bounded but unknown exogenous random disturbances. A 
smallest set refers to a smallest normed set of non-parametric uncertainties, i.e. , unmodeled 
dynamics. An illustrative example is given in section 5 which includes the computation and 
comparison of three model validating sets. Section 6 concludes the paper. 

2 Problem Definition 

For a given approximately linear and time-invariant physical system, suppose a controller, 
K , is known to be internally stabilizing and we can measure its inputs, u, command inputs, r, 
and the outputs, y , as shown in Figure 1. The exogenous disturbances consist of measurement 
noise at the output and exogenous disturbances either through the control input channel or 
through a separate path. To accommodate this, a noise filter, V, and a disturbance filter 
imbedded in the augmented nominal plant, P are assumed. Both filters are assumed driven 
by unknown but bounded independent random signals. The closed loop output of the plant 
model, y, is a sum of responses due to filtered noise, v, and disturbances, e, known input 
signal, r, all subject to feedback between a known controller, K , and a plant which belongs 
to a set defined by an upper LFT model, 

F U (P, A) := [P 22 , P 23 ] + P 2 iA(/ - Pi 1 A) _1 [P 12 , P 13 ] 
where the structured uncertainty is defined by A G P, where 

V := {A G C mxn : A = diag^, . . . , S r I nr , A r+1 , . . . , A T ), S t G P ( , A, f G C"‘ " } (1) 
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and r denotes the number of uncertainty blocks and T{ is either the field 1Z of real numbers 
or the field C of complex numbers at the designer’s choosing. A corresponding set of bounded 
structured uncertainty is defined by 


'Dw '■ — {A G H : A — A b W, ct(A^) < 1} 


( 2 ) 


where the elements of the diagonal matrix, 11' := diag( ic\ I ai . . . . , w t I Ut ), denotes the scaling- 
radii applied to the blocks of the structured uncertainty unit ball as defined by <r( A B ) < 1. 

For convenience, combine both exogenous disturbances e and v into a single disturbance 
vector 

* = { l } < 3 > 

The model in a canonical form is shown in Figure 1. For an open loop problem, G(P, K , V) = 
G(P, 0,1'), and r = u, and all the remaining developments remain the same. Of course 
the coupling effects of the noise and disturbances with the model uncertainties are more 
complicated due to feedback. Note also that since r and K are assumed known, as is 
typical in many applications, if the output error, e y := y — y = 0, then the input error, 
e u := u — u = 0. This means that it is sufficient to impose e y = 0 only. 

Suppose the measurements are taken in the discrete-time domain and consider a dis- 
crete frequency domain formulation, so called Constant Matrix case in [3]. For simplicity, 
we assume that a discrete Fourier transform has been performed and do not consider the 
additional affects of realistic signal conditioning operations typically performed on the raw 
discrete-time signals. Since <r(A) = max i (j(A i ) and we are primarily concerned with the 
size of the uncertainty blocks, it is of interest to note that uncertainty bound in terms of its 
maximum singular value can be written as a ratio of norms 


cr(A, ; ) := sup 


1 1 Ai^H ^ ||6 


VII — 


1 < i < t 


where col (6 .••••. 6 ) and col(//i, • • • , rj T ) are the partitioning of the vectors £ and rj which 
conform to the block diagonal partition of A in Equation (1). Of course fictitious signals £ and 


r) cannot be measured nor are they arbitrary so that it is necessary to look at their dependence 
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on real signals u and y and their transmission through a postulated system G(P, K, V) at 
each frequency. The signals, £ and r), whose norm ratios determine the uncertainty sizes, 
must be consistent with their transmission through the nominal plant model and uncertainty 
structure. The signals must also reproduce the measured inputs and outputs with some help 
from simulated exogenous disturbances which are assumed unknown but bounded by known 
bounds. 

The output error is given by 

(■y : !l // t", G 2 i£ — G 2 2 fi (5) 

where e° := y — G 23 r denotes the nominal output error due to the nominal plant in closed 
loop. The terms G 2 i£ and G 22 i3 in Equation (5) represent the uncertainty freedoms in an 
attempt to negate the nominal output error. The first term corresponds to the uncertainty 
freedom in an LFT structured uncertainty model, and the second term is an allowance given 
to unknown but bounded exogenous noise and disturbances in the system. 

Definition (Model Validation) [3]: 

Given measurements of the input signal, r, output signal, y, and noise filter, V, an 
augmented nominal plant model, P, with a disturbance filter, a controller, K , and a set of 
bounded structured uncertainty, V w . The set of plants (robust control design models) 

Vw '■= {F U (P, A), A G T>w} 


is said to be a model validating set if it contains an uncertainty model A G V w such that 
there exists exogenous disturbance signals, e, v with ||/2|| < 1 for which 


y = P u (G(PK 1 V),A){ v 

r 

In short, a set of plants is model validating if it can reproduce the given measurements 
while subject to a priori constraints, but of course, as noted earlier [3]-[6], one can never 
really “validate” a model since fresh data could potentially invalidate it. Notice that the 
assumption in the exogenous signal is different from the original problem definition in [3]. 
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A comment is in order about the combining of the two exogenous disturbance vectors 
e and v into the single vector 3. They are applied at dissimilar points of the plant model, 
and in many ways it would seem natural to keep them separate and, perhaps, define model 
validating set in terms of independent bounds on the two signals; typically, ||e|| < 1 and 
ll^ll < 1. However, the mathematics used in Lemma 2 in deriving the parameterization of 
model validating signals requires that ,3 be bounded in the ^ (Euclidean) norm, and the 
argument would not work if, instead, / 3 were bounded by the hybrid / ^ 2 ) norm which 
would be equivalent to the individual bounds on e and u. However, the two conditions are 
fairly close to each other: 

1 1 3 1 1 < 1 => ||e|| < 1 and ||zz|| < 1 => ||/?|| < ypi 

In the case that one or the other of e and v is negligible, then ||/j|| is approximately equal 
to the norm of the other, and even in the general case, the discrepancy between the two 
bounding paradigms is only y/2, independent of the dimensions of v and e. 

In this paper, we characterize model validation in a way that will allow a convenient pa- 
rameterization of all model validating uncertainty sets with respect to available input/output 
measurement for a given LFT structure. The idea is that the engineer will ultimately select 
a particulat set of model validating uncertainty weights based on the particular application. 
To this end, we first investigate a feasibility condition (or necessary condition) for model 
validation: at each frequency, is there a pair (£,/?) subject to ||/3|| < 1 which makes the out- 
put error in Equation (5) zero? Note that once G, r, 3. and f are specified, rj is completely 
determined. The next step is to incorporate the constraints due to a priori structure in the 
uncertainties which may limit the feasible (£, rj) signals. This leads to necessary and suffi- 
cient conditions for model validation. The final step is to parameterize all model validating 
sets of plant models. 

It is clear that if the noise vector v is not restricted by a fixed bound, then any output 
residual can be zeroed out (without any help from A and/or e) if the noise filter V is 
non-singular. In this paper we assume that the noise filter V is given (as part of the a 
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priori model assumption or a reasonable model determined from earlier system identification 
experiments) and the noise vector at each frequency is norm bounded by 1. This output 
noise model can be viewed as a model of a broad band exogenous noise typified by sensor 
noise. Of course when judicisoulv chosen, V can reflect a priori bounds on the noise intensity 
or power spectrum of the unknown exogenous signal over a bandwidth of interest . The input 
vector, e, is meant to denote unknown but bounded exogenous disturbances on the physical 
system. As an example, in aeroservoelasticity applications, e will denote white noise input 
to a gust filter imbedded in P\ 2 and P 22 . 


3 Parameterization of Model Validating Sets 


In this section, we develop a theory to effectively parameterize all model validating uncer- 
tainty sets which satisfy a priori assumptions on the LFT structure and exogenous distur- 
bances. We begin by first addressing the question: when does there exist a norm bounded 
combined exogenous disturbance (3 with ||/i|| < 1 and a £ £ C n z such that e y = 0? At this 
point, notice that £ is not required to be limited by any given bound or structure. To answer 
the above question we first state a lemma, set some notation, and make some observations. 


Lemma 1: 

Let A be a matrix whose singular value decomposition (SVD) is 


A = [ Ui U 2 ' 


Si 0 
0 0 


H 


V{ 

v? 


( 6 ) 


where U = [Ui U 2 \ and V = [Iq V 2 ] are unitary, is diagonal and nonsingular, and the block 
matrix partitionings are conformal. Let b be a vector and c a non-negative real constant. 
Then the inequality 


has a solution x if and only if 


| . hr + b\\ < c 


\\U?b\\ < c. 


( 7 ) 

( 8 ) 
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Then the general solution to the inequality (7) is parameterized by 


( 9 ) 

( 10 ) 


• = Vi y + l 
y = E ^(w-Ufb) 

where z is arbitrary, and w is any vector with 

\H\ < y/<? ~ \\U?b\\*. (11) 


Remark: 

One or both of two special circumstances may apply to the SVD of A: A might be of 
full row [column] rank in which case the U 2 [17] matrix is missing from the SVD, as is the 
bottom row [right column] of zeros in the diagonal matrix of singular values. The following 
modifications make Lemma 1 and its proof remain correct under any combination of these 
circumstances. If A is of full row rank, then consider ||L r | J 6|| to be 0 and remove the U 2 lines 
from Equations (12) and (13). If A is of full column rank, then remove any reference to ,2 
and 13 and realize that w subject to the inequality (11) is the only free parameter in the 
parameterization of solutions to the inequality (7). The necessary and sufficient condition in 
equation (8) is saying that the norm of the projection error of b onto lh must be less than 
or equal to the given inequality bound. 


Proof of Lemma 1: 

Note that x can always be written as \\y + V 2 z by taking y = Vf 1 x and z = Vf 1 x. Now, 
from the SVD of 4, AV 2 z = 0, and A Vj y = U\'h\y. Therefore, for arbitrary y and if 
x = V] y + V 2 z then ||4a; + b\\ 2 < c 2 if and only if \\Ui'Eiy + 6|| 2 < c 2 . Since U is unitary, 
multiplying by U H preserves norm, so the last inequality is equivalent to 


H 


U{ 

u 2 


H 


In its turn, this is equivalent to 


U\P\y + 


Si y + U?b 
U»b 


U{ 

U? 


< c 2 . 


< c 2 . 


( 12 ) 


(13) 
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This can be rewritten as 


\\X iy + U*b\\ 2 <(?-\\U*b\\ 2 . (14) 

So on the one hand, if inequality (7) has a solution, then that solution fits the param- 
eterization given in Equations (9) - (11) and condition (8) is satisfied. On the other hand, 
if condition (8) is satisfied and first w and z and then y and finally x are picked according 
to the parameterization given in Equations (9) - (11), then this x provides a solution to 
inequality (7). □ 

3.1 Constant Matrix Test 

If we set M := [G 2 i,G 22 ], then the condition e y = 0 can be written as 

This provides the first necessary condition for the existence of j3 with ||/3|| < 1 and £ for 
which e y = 0; namely: 

e° e Im(M) (16) 

Observe that if V is invertible, then condition (16) is true, since M has full row rank. 

Physically, this means that if the noise model at output is allowed to influence all output 
channels, then any output signal can be validated (without any help from A and/or e) if the 
noise is not constrained by a bound. Also, condition (16) is necessary and sufficient for the 
existence of [3 G C ne+n v and £ for which e y = 0. The remainder of this discussion is aimed 
at finding a condition to insure that ||/3|| < 1. 

If condition (16) is satisfied, then Equation (15) is solvable, and a complete parameteri- 
zation of the solutions is given by 

| | | = M + el + N m 9, (17) 

where N M is a matrix whose columns form a basis for Ker(M), and the parameter 0 is 
arbitrary. The notation (•)+ denotes the Moore-Penrose pseudo-inverse of (•). 
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Notation: Let a be a vector defined by an equation of the form a = Ax + By + ■ ■ and let 
b be a subvector of a. Then by {A)b, ( B)b , etc., we are referring to those submatrices of A, B , 
etc., containing the complete rows of the respective matrices so that b = ( A)i,x + {B)t,y + • • •. 
Equation (17) then implies that 

fi=(M + ) p el + (N M ) p e (18) 

and the problem is to characterize those 0 for which ||/3|| < 1. 

This can be done by a direct application of Lemma 1. First perform an SVD of ( N M )p 
as in Lemma 1: 

(w,a = [ a n ] [ E 0 > l ] [ ] as) 

Then, by Lemma 1, Equation (18) can be solved for ft with ||/3|| < 1 if and only if 

\\T?(M*) e e;\\ < 1. ( 20 ) 

Then, the general solution to ||/3|| < 1 subject to Equation (18) is given by 

6 = l']'- I l -2‘P (21) 

7 = (22) 

where ft G C n,:: is arbitrary, and ft G C n ° is any vector with 

iqil < K := v'l - . (23) 

In light of the definition in Equation (23), the necessary condition in (20) implies b 0 > 0. 
Lemma 2 (Constant Matrix Test): 

With the context and notation established in the previous paragraphs, the following are 
equivalent: 

1. There exists {ft ft), ||/3|| < 1 such that y = G 2 i£, + G 22 ft + GW’- 

2. Conditions (16) and (20) hold. 
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If these conditions are satisfied, then all such triples (£, e, v) are given by the parameterization 


£ 

P 


c 0 

fio 


n 


i' 


where 


SI 

0 


-1 


O := N m U 

:= [A/' .V A/ ((.\ A , ),)•(.!/•); <‘; 

and where r W is arbitrary, and (p satisfies inequality (23). 


Co 

fio 


(24) 


Remark: The number of free parameters (</>, V’) as defined by Equation (24) is given 
by n f /) + n.0 = dim(Ker(M )), and the dimension of the identity matrix, /o. is n 4) + — 

rank(N M ) g . 


Proof of Lemma 2: 

The existence of a triple (£, e, y) with the desired properties is equivalent to the solvability 
of Equation (15) with ||/3|| < 1. We have already seen that Equation (15) can be solved with 
||/?|| < 1 if and only if conditions (16) and (20) hold, and all such solutions are parameterized 
by Equations (21) - (23). When a value of 9 as given by the parameterization in Equations 
(21) - (23) is substituted in Equation (17) to eliminate 9 and then 7, the expression for the 
triple (£, e, u) given Equation (24) results. □ 


Lemma 2 gives a test resulting in either a yes or no answer. It is only concerned with 
testing the richness of the a priori LFT uncertainty structure and chosen levels of measure- 
ment noise and disturbance allowance against a given set of measured input and output data. 
Whether such f can be generated through the LFT uncertainty and how large it must be 
remains to be seen. If the test in Lemma 2 fails, then the model is invalidated either due to 
overly restricted levels of noise and/or disturbance and/or insufficiently rich uncertainty LFT 
structure. What course of action should be taken if the Lemma 2 test fails is not considered 
in this paper. Of course, at this point, increasing the magnitude of the LFT uncertainty 
bound will not help. 
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Suppose the conditions of Lemma 2 are satisfied. Indeed, Lemma 2 gives a parameteri- 
zation of the set of all £ that produces zero output error. That £ is given by this parameter- 
ization provides a necessary condition that £ be a signal in a model validated robust control 
design model. For sufficiency, £ must also satisfy the G — A feedback conditions 


A 7) 

(25) 

Gn£ + Gup + Gi^r 

(26) 


Since 7) can be readily computed from Equation (26) for a given £ and /!. we group Equations 
(24) and (26) as follows: 


V 

£ 


£ 0 


Gfii G 1 2 

Ori£X(n e +n„) 



(27) 


Here, 

Vo '■= [ Gfil <-712 ] | ^° | + ^13^1 

and the norm of (j) is subject to condition (23). Equation (27) characterizes the set of all 
(£, 7 j) vectors that produces zero output error. Of course, this set may be further constrained 
by the uncertainty structure given by Equation (25), which motivates the next lemma. 
Consider a basic fact from linear algebra as noted earlier in [8]: 

Lemma 3: 

If u E C m \ v E C n , v ^ 0, then there exists A E C mxn such that Av = u, and a(A) = jj^jj. 

Remark: If Av = u, then ||H|| = a (A) > ||^||, so this lemma is asserting that an A of 
the minimal possible norm does exist which maps v onto u. This will find application in 
this paper in demonstrating the existence of model validating blocks A,; of minimal possible 
norm. 


3.2 Full Complex Blocks Only 

We now work toward determining when there exist model validating A E V with the implicit 
assumption of structured, full complex blocks only for A in this section. To this end, we 
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start with uncertainty signals £ and rj which satisfy Equation (27) and look for A for which 
Equation (25) is also satisfied. Partition £ and r) into components corresponding to the block 
structure of A so that Equation (25) could be written as 

& A,//,. * = (28) 

Since, for each i = 1 , . . . , r, £* = A ; // ; . one can never have both ^ 0 and ^ = 0 (cf. 

Lemma 3). However, the parameterization in Equation (27) does not guarantee that these 
conditions never occur. To have some terminology to use to indicate that we are excluding 
this possibility, we make the following definition: 

Definition (D -realizable): 

A signal pair (£, 77 ) will be called ^-realizable if, for each i = 1, . . . , r, either = 0 or 
Vi ± 0 . 


Satisfaction of the test in Lemma 2 allows a cancellation of the nominal output error 
by a combination of the fictitious signals from the uncertainty block and a norm bounded 
exogenous disturbance while it will be shown that Lemma 3 guarantees that a structured, full 
complex uncertainty always exists for any D-realizable pair (£, r\) which satisfy Equation (27). 
We state an existence condition and a parameterization of all model validating uncertainty 
sets as follows: 


Theorem 1 (structured, full complex blocks): 

Suppose the conditions given in Lemma 2 are satisfied. Then, all model validating sets 
of plants are given by 

Vwh, '■ {F U {P, A), A G V w } (29) 


where fi G C n * , (f) G C n ° . ||0|| < b 0 , W := diag(uq/ ni , . . . , w t I U t ) is any diagonal complex 
matrix whose diagonal elements satisfy 


Wi > 



7 ' = 1, . . . , r, 


(30) 
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and the (£, 77 ) pair parameterized by <p and ip as given in Equation (27) is ©-realizable. 
Proof of Theorem 1: 

It is first demonstrated that each Vw^tp described in the statement of Theorem 1 is a 
model validating set. Since the conditions of Lemma 2 are satisfied, the f given by Equation 
(27) combines with an exogenous disturbance f3 for which ||/3|| < 1 to produce e y = 0. 
Therefore, Tw^ will be model validating if there exists A e V w for which Equation (28) 
is satisfied. We construct this A block by block. If = 0, then A< and A Bi may be taken 
to be 0 and a(A Bi ) = 0. If ^ 0, then since (£, 77 ) is a ©-realizable pair, f], ^ 0 and by 
Lemma 3, there exists a A, ; with A * 77 ,; = and <r(Aj) = ||£j||/|| 77 j||. This means that if 
Wi 7 ^ 0 and A Bi := (l/iy^Aj, then a(A Bi ) < 1 (if w, = 0, A Bi may be chosen arbitrarily 
with a{ A Bi ) < 1). It follows that if A B := diag(A B1 , A B2 , • • • , A Sr ), then a(A B ) < 1 . Then 
A = A b W £ ©iv . so ©iv vc is model validating. 

Now let ©iv be an arbitrary model validating set. Let A £ V w zero the output error for 
some J3 where ||/3|| < 1. Then there exist signals f and 77 which satisfy Equations (25) and 
(26) such that f zeros out the error with some d where ||/3|| < 1. Then by Lemma 2, there 
exist (j) satisfying condition (23) and r tp such that f and i] are expressed in terms of (p and r tp 
by Equation (27). Since = A,;/, for all i = 1, . . . , r, the pair (£, 77 ) must be ©-realizable, 
and ||£j|| = || A* 77 * || < ( A* ) 1 1 77 * 1 1 < \wi \ 1 1 77 * 1 1 , the last inequality holding because A e V w . 

This establishes that Vw = Vwh>- n 

For practical applications where the ©-realizable assumption is typically satisfied, Lemma 
2 can be viewed as a necessary and sufficient condition for the existence of a model validat- 
ing set for an LFT with only structured full complex blocks. The issue that remains is how 
large must the uncertainty size be for model validation, which is addressed in Theorem 1. 
For fixed parameters <p and ip , the smallest model validating uncertainty set is given by 
w t = j|| L || . 1 = 1, . . . , r. This is the minimum norm model validating result reported in [ 8 ] 

and any choice of uncertainty set bounds { i = 1 r) for which I wA < inf j|^4| for 

m\»o Hwi 

some t. will not be a model validating set. 
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3.3 Repeated Scalar Blocks 

For a more general uncertainty structure which includes repeated and/or real scalar blocks, 
we assume for convenience, that all repeated scalar blocks are grouped into the first r blocks 
in A as given in Equation (1). Since A is further restricted, the conditions in Lemma 2 
or Theorem 1 becomes only a necessary condition for model validation. So, with repeated 
scalar blocks, we ask whether, among all D-realizable pairs (£, rf) satisfying Equation (27) 
with <f> subject to the norm condition (23), a pair exists for which a A of the form given in 
Equation (1) also exists so that Equation (25) is satisfied? 

First, note that given any such D-realizable pair (^, 77) , A,, +1 , . . ., A r always exist by 
Lemma 3, so we need to consider only the existence of the first r blocks, S 1 /„, , . . S r I nr where 
Sj, e Tj. Let us denote n := Yi-i n i = ™ '■= YJi \ m i- A model validating uncertainty set 
exists for a system with repeated scalar block if and only if there exists S, G JF) , i = 1, . . . , r, 
'A </>, ||</>|| < b 0 such that 

A Sj i/;. i = 1 , . . . , r (31) 

where from Equation (27): 

f. = fo.. + a { / j (32) 

ij = + [Cn O j] , - - 1 / 1 (33) 

The subscript i indicates that the correct blocks of rows have been selected for Equations 
(32) and (33) to make sense in the context of Equation (27) and the decompositions of £ 
and 7]. 

The condition in Equation (31) can be seen as a collinearitv condition in the vector 
space C n ‘ with coefficients from the field JF). Consequently, a measure of distance between 
two subspaces can be used (see for example [12]) 

dist^ ) K„l ) .):=||Pf‘ l _pW)|| (34) 

where and Pjjf^ denote orthogonal projections onto the subspaces spanned over the 
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field T{ by the single vectors £* and rji, respectively. As an illustration, in the case of one 
dimensional subspace in 7Z 3 , the distance measure as given by Equation (34) reduces to the 
sine of the angle between the two lines. We summarize our results as follows: 

Theorem 2 (with repeated scalar block): 

(a) Suppose the conditions in Lemma 2 are satisfied. Then a model validating set exists 
with A i = 5jl nn i = 1, . . . , r, 6, G T\, if and only if there exists A, and (]) with ||0|| < b 0 such 
that the (£, rj ) pair as given in Equation (27) is D-realizable and for each i = 1, . . . , r, 

either ^ = 0 or dist^(£*, r/i) = 0 (35) 

where (, and /],; are given by Equations (32) and (33). 

(b) Furthermore, if a model validating set exists, then all such sets are given by Equation 
(29) where W is any diagonal complex matrix whose diagonal elements satisfy condition (30) 
and the (£, 77 ) pair parameterized by and r tp as given in Equation (27) is ^-realizable and 
satisfies condition (35). 

In order to prove Theorem 2, we first introduce a lemma. 

Lemma f: 

For each fixed %. condition (35) holds if and only if there exists a Si G T, such that 
h,i //; . 

Proof of Theorem 2: 

First, suppose that there exist 'W and <p with ||<(>|| < b 0 such that the (^, 77 ) pair as given 
in Equation (27) is P-realizable and condition (35) holds. Then by Theorem 1, any set 
Vwtfiip as given in Equation (29) with W satisfying condition (30) is model validating. Let 
A G V w be a model validating uncertainty model. The key properties of A are that, for 
each i = 1 ,...,r, ||A,;|| < |iuj| and fi = A,r/j. Because of condition (35), Lemma 4 tells 
us that we can replace each A* for i = 1 by a matrix of the form 4 ;J n . with <5* G 
and still have & = A,;^. Also, since & = 5^ and W satisfies condition (30), we still have 
IIAiH < | Wj | with this new A i: , This establishes the “if” part of Theorem 2(a), and shows 
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that every model validating set with the specified repeated scalar blocks has the form shown 
in Theorem 2(b). 

Now suppose that a model validating set V exists which contains a model validating- 
uncertainty model A with Aj = i = 1, . . . , r, d, : E Tj. Theorem 1 tells us that V must 

satisfy all of the conclusions of Theorem 2(b) except possibly the condition (35). However, 
since for each i = 1, . . . , r, A, = <5* it follows from Lemma 4 that condition (35) is satisfied. 
This completes the proof that every model validating set falls under the description given in 
Theorem 2(b), and completes the “only if” part of the proof of Theorem 2(a). □ 

4 Uncertainty Bound Tradeoff 

With the parameterization of all model validating uncertainty sets given by Theorems 1 
and 2, a controls engineer still faces the issue of how to utilize the remaining freedom. 
Specifically, one may ask: can we find a smallest set from the given parameterization of 
all model validating uncertainty sets? This begs the issue of which “smallest set” for this 
underdetermined problem. For a single uncertainty block problem, a smallest-norm model 
validating uncertainty appears to be a physically reasonable uncertainty bound based on 
Ockham’s razor argument in modeling physical systems. However, for problems with a 
general LFT uncertainty structure which can be viewed as a multi-objective problem, an 
uncertainty bound having a smallest norm may not have any concrete physical justification. 
This is because for problems with multiple uncertainty blocks, their relative numerical values 
do not necessarily indicate their relative physical significance. For example, in robust stability 
[13], the determination of whether a controller guaranteeing robust stability exists or not 
may depend more strongly on the distribution of the uncertainty bounds over a given set of 
uncertainty components than on the size of largest uncertainty component. 

In this section, we outline two algorithms based on constrained nonlinear optimization 
to determine a smallest uncertainty set. One starting point for an optimization would be 
for the designer to select an augmented nominal plant model, P, a noise allowance filter, V. 
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and a matrix, W, of desired levels of uncertainty in the blocks. The important feature of 
W is that the relative sizes of the Wj, reflect the designer’s wishes as to the relative size of 
the uncertainty levels in the different blocks. A smallest uncertainty set search algorithm is 
outlined as follows: 

1. select (P,V,V) 

2. test feasibility of (P,V,V) against (r, y). 

3. select W 

4. find smallest x such that V( x w) is model validating 

If the feasibility test in Step 2 fails, redo Step 1 until it passes. For uncertainty structure 
with only full complex blocks, Step 2 is simply a constant matrix test given in Section 3.1, 
otherwise, Step 2 can be included as a constraint in Step 4 (see Remark in Section 4.3). Using 
V(xw) from Step 4 and a chosen performance weight, a robust controller can be obtained 
by n synthesis. If the closed loop system is determined to be unsatisfactory or suspected to 
be unrealistic, a new model validating set can be obtained by reselecting the weights W in 
Step 3 and repeating Step 4. In Step 4, if the chosen W is a desired set of weights and the 
computed smallest x is less than 1, then the postulated robust control design model, V w 
can be considered to be model validating. 

As a prelude to trying to optimize uncertainty levels during a tradeoff, the necessary 
conditions (16) and (20) for the existence of model validating sets, should be checked. If all 
of the output channels are being modeled as having noise in them, so that the V matrix is 
non-singular, then the matrix M is of full rank, and condition (16) always holds. However, 
if the diagonal matrix V does have zeros on the diagonal (or, for the general case V > 0 
for fully populated V), then a check on condition (16) can be made by first performing an 
SVD on M: M = l u l '// where U M and V M are unitary and E M is a non-negative real 
diagonal matrix of the same shape as M whose diagonal elements are in decreasing order. 
Condition (16) also holds if M is full rank which is equivalent to S M having no zero rows. 


19 



If M is rank deficient, then partition U M = [Umi,Um2\ where the block U M i corresponds to 
the non-zero rows of S M and the block Um 2 corresponds to the zero rows of S M . Then a 
necessary and sufficient condition in the case of rank deficient M for condition (16) to hold 
is that U M 2e° y = 0. To verify condition (20), the SVD of M is computed in the form: 


M — ( d/i Um 2 


N - \i i 0 
0 0 


V H 

v Ml 

v H 

1 M2 


(36) 


In this decomposition, U = \Umi Um2\ and V = [Vmi Vm 2 ] are unitary, S M i is diagonal 
and nonsingular, and the block matrix partitionings are conformal. Then, in Equation (17), 


M + = {/i and N M can be taken to be V M 2- Then, (N M )p is decomposed as 


in 


Equation (19), and all of the components are at hand to execute the test in condition (20). 


4.1 Full Complex Blocks 


The idea is to normalize all uncertainties using desired levels of uncertainty in the uncertainty 
blocks, and then seek the smallest model validating scaled set. Specifically, we propose using- 
nonlinear constrained optimization with ip and (p as the design parameters to find a minimal 
positive x such that V( x w) is a model validating set. 

By Theorem 1, V( x w) is a model validating set if there exist ip and <p with ||0|| < b 0 such 


that 



(37) 


where the vectors £ and ij calculated from ip and (p in Equation (27) form a 'P- realizable 
pair. This implies that V( x w) = 'P(xW)tM)- By squaring and clearing fractions, the previous 
inequality can be combined with the 27-realizability condition in the single inequality in 
Equation (39). This also has the advantage of being a polynomial in x and the components 
of p and rj. The optimization problem can now be stated: 


A Smallest Set (Full Complex Blocks): 


mime (38) 
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subject to 


£ry 

to 

1 

H 

to 

jT 

to 

to 

< 

0 , / 1 ... 

• ,T 

(39) 

X 

> 

0 


(40) 

M 

< 

b 0 


(41) 


where / and //, are given by Equations (32) and (33). 

For the special case where the noise is known or given, the parameter <p is unbounded 
and Equation (41) is no longer needed. This is the case derived earlier as a minimum norm 
model validating solution [8]. For the class of problems where the uncertainty structure is 
such that all r/, 's are completely determined by input and output data (for example additive, 
multiplicative, additive feedback, etc), a feasible set in the above optimization algorithm 
is convex for any given x. Hence for this special class of uncertainty structure, the model 
validation test (with x = 1) can be seen as a convex feasibility problem as noted earlier [4], 

4.2 Repeated Scalar Blocks 

For the case with repeated and/or real scalar blocks, an optimization algorithm similar to 
full complex block case but with the additional collinearity condition is proposed. Similarly, 
by Theorem 2, V( x w) is a model validating set if there exist ip and <p with ||0|| < b 0 such that 
condition (35) is satisfied where and i], is a D-realizable pair parameterized by Equations 
(32) and (33). Instead of using the distance condition in Equation (35) to guarantee existence, 
Equation (31) is used. The tradeoff is that Equation (35) leads to a quartic in the design 
variables while Equation (31) leads to only a quadratic at the expense of additional variables, 
hi,. . ., S r . Note that the collinearity condition for the set of r repeated scalar blocks and the 
D-realizability condition leads to a simplification of the first r set of inequality constraints 
in Equation (39). The optimization problem can now be stated: 

A Smallest Set (With Repeated and/or Real Scalar Blocks): 
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min x (42) 

i>,<p,Si„.„5 r ,X 

subject to 

Si e Ti i 1 (43) 

6 = i = (44) 

\5i\ 2 - x 2 \Wi\ 2 < 0, i 1 ,r (45) 

||6|| 2 -x 2 K| 2 ||%|| 2 < 0, i = r + l r (46) 

x > 0 (47) 

M < b 0 (48) 

where & and //, are given by Equations (32) and (33). 


Remark: In Section 4 constant matrix tests were given which were necessary, but not 
sufficient for the existence of a model validating set for the case that some uncertainty 
blocks have repeated scalar blocks. Such a sufficient condition is found by the existence of a 
feasible point in the preceding optimization problem. In particular, in order to have a model 
validating set in this case, it is necessary to satisfy constraints (43) and (44). The existence 
of a feasible point is also sufficient for satisfying ©-realizability. Note also that the above 
constraints involve polynominals in </>, ■</>, di, . . . , S r , x which are at most quart ic. 

To summarize, the above optimization algorithm has the following physical significances. 
The cost in Equation (42) represents a positive scaling factor of the normalized (by user 
provided desirable weights) uncertainty norm bounds for each component. Inequalities (45) 
and (46) represent the scaled bounds on r repeated scalar uncertainties and the r — r non- 
repeated full complex uncertainty bounds, respectively. A violation of these inequalities 
implies that x\wi\ is not an upper bound on the ratio of signal norms, i.e., it fails as an 
uncertainty bound. Clearly, inequalities (45) and (46) will more likely be satisfied with 
larger weights x\wi\, which makes intuitive sense. The collinearity condition in Equation 
(44) represents the necessary structural constraints due to the repeated scalar uncertainties. 
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Inequality (47) is the non-negative condition on the uncertainty scale factor. Note that 
x > 1 indicates that the current scaling makes the uncertainty bounds larger or equal to the 
a priori target while x < 1 indicates that there exists a smaller (for every component ) model 
validating set than the a priori target. Finally, inequality (48) represents a noise allowance. 

4.2.1 Real Parametric Uncertainties 

In the uncertainty optimization algorithm such as Equations (42)-(48), optimistic (to an 
unknown degree) levels of non-parametric uncertainty bounds are expected since the uncer- 
tain parameters are free to be a function of frequency although subjected to fixed lower and 
upper bounds. In ^ analysis and synthesis, the uncertainty is only restricted to be norm 
bounded and stable and therefore will lead to conservative results. The net result is that 
an optimistic prediction of non-parametric uncertainty bounds will be compensated to some 
unknown degree by the conservative prediction from /u analysis and synthesis. 

Suppose the repeated and/or real scalar blocks are independent of frequencies, i.e. , the 
parameter is unknown but constant and we wish to accomodate this by specifying lower 
and upper bounds to design for anticipated changes. Of course parameter estimation can 
precede the specification of the above bounds. Instead of jointly optimizing a set of uncer- 
tain parameters with non-parametric uncertainties as shown in Equations (45) and (46), we 
propose solving for the smallest non-parametric uncertainties subject to a priori allowance 
in exogenous disturbance/noise and parametric uncertainties. This can be implemented by 
eliminating the scale factor x in Equation (45). A practical justification for this approach is 
that (physical) parameters are generally better modeled then non-parametric uncertainties 
since the latter type typically originates from complex dynamics with many details which 
are difficult and subsequently not explicitly modeled, viz “unmodeled dynamics”. There- 
fore, to determine a smallest model validating uncertainty set for problems with parametric 
uncertainties, we propose (a) introducing sufficient non-parametric uncertainties to satisfy 
model validation conditions (constant matrix tests) and then (b) introducing varying levels 
of parametric uncertainty allowances to tradeoff with smallest non-parametric uncertainty 
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levels. 


5 Illustrative Example 

5.1 Simulated System and Nominal Model 

Consider a simulated (true) system composed of two lumped masses-spring-dashpot with 2 
inputs and 2 outputs which are fixed at two ends as given in [7]. For the purpose of this 
study, suppose the nominal plant model is chosen to consist of only the first mode to reflect 
a reduced order model. The neglected second structural mode reflects unmodeled dynam- 
ics. The nominal plant model consisting of the single mode is assumed to have parametric 
uncertainties in the frequency and damping. The two structural resonances corresponds to 
approximately 2.9 and 5.0 Hertz with damping ratios of 6 and 10 percents respectively. 

Figure 2 shows the Bode plot of the true system, nominal plant model, and their differ- 
ence. Notice that the largest differences occur at the structural resonances. The first peak 
in the difference is due to the parametric uncertainty in the first mode whereas the second 
peak is due to the unmodeled dynamics corresponding to the truncated second mode. 

With the chosen nominal plant and the uncertainty structure, the system interconnection 
is constructed. The frequency and damping parametric uncertainties are modeled in discrete 
time state-space while the non-parametric additive uncertainty is between the control inputs 
and plant outputs. The frequency and damping parametric uncertainty in the first mode as 
represented in the discrete time domain leads to a repeated parametric uncertainty of order 
2. A set of 1024 time points are assumed sampled at 50 Hz. The simulated measurement 
noise and input test signals are generated by filtering independent uniform random signals 
through two different wideband filters. The simulated output signal to noise ratio ranged 
from 1 to 10 3 over bandwidth of interest. 
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5.2 Parametric Study 


Consider the problem of finding the smallest unmodeled dynamics uncertainty for model 
validation subject to bounded noise and parametric uncertainty. A smallest set algorithm 
defined by Equations (42) - (48) is applied with a modification to Equation (45) as discussed 
in Section 4.2.1. Analytical sensitivity formulae were used with a constrained optimization 
routine in [14] based on a quasi-Newton algorithm. Since unlimited additive uncertainty is 
included, existence conditions given in Lemma 2 are satisfied. We consider nine cases (see 
Table 1 ) corresponding to a combination of 3 levels of parametric uncertainty allowance and 
3 levels of output noise allowance. In the determination of a smallest model validating set, 
parametric uncertainties are allowed to be functions of frequency as long as they satisfy a 
constant norm bound over frequency. As discussed earlier, this will lead to an optimistic level 
of non-parametric model validating uncertainty. However, in the actual design of controllers, 
the a priori parametric norm bound is used for the robust control design model, instead of 
the smaller converged parametric uncertainty value. 

5.2.1 Output Noise Allowance 

As shown in Table 1, three levels of the output noise allowance are considered. The peak 
2-norm of the true noise spectrum over the bandwidth of interest, 1 to 10 Hz, is denoted by 
Vpe ak ■= 0.078/2x2- In Cases la, 2a, and 3a, a high noise allowance corresponding to V pea k is 
assumed, while a moderate noise level of V rms = diag(.0321, .0332) is assumed in Cases lb, 
2b, and 3b. This root-mean-square over the above bandwidth is computed independently 
for each output channel. Finally, a low noise allowance of 1% of V pea k is assumed in Cases 
lc, 2c, and 3c. In the actual model validation computations, since the allowed independent 
white noise, u, is bounded by unity, the noise filter is scaled by a factor of \/2. 

5.2.2 Eigenvalue Uncertainty Allowance 

In Cases l(a-c), an optimistically small level of parametric uncertainty allowance is assumed, 
Cases 2(a-c) represents a realistic or ball park estimate, and finally Cases 3(a-c) represents 
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overly pessimistic levels of the actual parametric uncertainty levels. Given parametric un- 
certainty bounds 

= -913 A .348?' + <5i + (^2?, |(5i| < hi | <^2 1 A S 2 

where 

X true = (.918 ± .358?, .756 ± .550?) 

find smallest model validating additive uncertainty. In other words, given specific output 
noise allowance, to what extent does specified parametric uncertainty levels contribute in 
reducing the additive errors? 

5.3 Smallest Model Validating Non-parametric Uncertainty 

Figure 3 shows the smallest levels of unmodeled additive dynamics needed for model val- 
idation for corresponding levels of parametric uncertainty and noise allowances with the 
following observations: 

• For Cases 1 (a-c) where a small level of parametric uncertainty allowance is given, 
the minimum model validating additive uncertainty levels show two peaks at the reso- 
nances. A large level of additive uncertainty is required to cover the parametric error 
in the first mode. 

• In Cases 2 (a-c), a ball park level of parametric uncertainty allowance helps to eliminate 
the additive uncertainty peak in the neighborhood of first mode. 

• An increase in the parametric uncertainty allowance as in Cases 3 (a-c) further reduces 
the additive uncertainty in the neighborhood of first mode. However, the minimum 
levels of additive uncertainty for model validation around the second mode frequencies 
did not drop significantly in spite of a large parametric uncertainty allowance. This is 
expected from a structural dynamics viewpoint whereby allowing changes in a modal 
parameter will not significantly affect other well separated structural resonances. 
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• Increasing noise allowance (up to the true peak noise value) has a negligible effect on 
the level of minimum model validating additive uncertainty although a slight drop in 
the uncertainty levels are noted. A noise allowance level of at least 100 x the peak noise 
was necessary to cover the entire output nominal error due to parametric uncertainty 
and unmodeled dynamics. 

5.4 Robust Performance Comparison 

Stable, low-order (up to 4th order) rational functions are used to over bound the model 
validating additive uncertainties shown in Figure 3 for the three cases using a routine from [2], 
Based on the three sets of fitted additive and corresponding parametric uncertainties, D-K 
Iterations with constant D-scales were used to design suboptimal // and skewed-// controllers, 
K 2 , and K* . and Kj Js , K 2 , and K* s , corresponding to Cases 1 to 3. The skewed-// 
controllers were obtained by computing a scaled sequence of // synthesis problems (see for 
example, [1], Chapter 8). The performance weight chosen is, W per j = 2/ 2X 2- All calculated 
lower and upper bounds were close which means that the bounds are a good estimate of 
the actual values. The skewed-// controller represents the optimal worst case performance 
controller over each model validating set. Although the uncertainty parameters are assumed 
and computed as real numbers, for simplicity they were assumed to be complex perturbations 
in the controller synthesis so that it is biased (larger //) towards Case 2 and more against 
Case 3. The performance indices for the controllers are summarized in Table 2. Case 3 
appears best (peak //, peak skewed-//), despite the assumption of d, G C. All three ii r , ea k 
occured around the first structural mode. 

As a final comparison, the true closed loop response is computed for all six controllers. 
The maximum singular value frequency responses of the simulated (true) closed loop transfer 
function matrices are shown in Figure 4. All three // controllers show equal performance 
on true model whereas Skewed-// controllers performed differently ( K 2 _ . K^_ show equal 
performance while gave poorer performance) Interestingly, this better performance (in 
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terms of //. Skewed-^, True response) occurs in spite of the fact that the uncertainty set 
used in K" 1 and were more conservative since parametric uncertainties were treated as 
complex in the controller synthesis and K" 1 and assumes larger parametric uncertainties. 

In summary, non-physical model validating sets (Cases 1,3) can produce as good (Case 
3) or worse result (Case 1) than a physically correct set (Case 2). In this particular example, 
there appears to be a preference of parametric uncertainty over additive uncertainty for 
better robust performance. This may be related to the fact that the open loop response 
shows a significantly larger first mode response than the second mode (see Figure 4) and 
the fact that the first mode is limited by parametric uncertainty which Case 1 practically 
ignores. 

6 Conclusions 

For models of physical systems where the uncertainty is described by a linear fractional trans- 
formation with unknown but bounded exogenous disturbances, conditions for the existence 
of a model validating set are derived. For the case with only structured full complex blocks, 
this condition is necessary and sufficient and can be readily tested. This test is actually a 
test on the richness of the a priori uncertainty structure enhanced by a frequency weighted 
unknown but bounded exogenous disturbances and measurement noise allowance. It is sig- 
nificant that this test applies for an arbitrary number of full complex blocks. For the more 
general case when repeated and/or real scalar uncertainties are also present, we have shown 
that an additional condition involving a collinearity test is required. Model validating sets 
of uncertainties are in general highly nonunique so that a useful design tool should provide 
an efficient way of trading off parametric, non-parametric and exogenous noise uncertainty 
levels. Based on a parametrization of all model validating sets, an optimization algorithm is 
proposed which seeks a smallest model validating set of non-parametric uncertainties while 
subject to fixed levels of parametric uncertainties and noise allowances. 
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Table 1: Parametric uncertainty and noise allowances. 



Si 

— 

1/ "s/2 * 1 noise 
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.01 


Case la 

.001 

.001 

^ peak 

Case lb 
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.001 
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•01 * 1 peak 
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Table 2: Peak values for //, and Skewed -/i for robust performance. 
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Figure 1: Uncertainty bound identification framework (top), in canonical form (bottom). 
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Figure 2: Bode plot of true system (solid), nominal model (dash), and difference (dot) for 
input 1/output 1. 
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Figure 4: Simulated (true) closed loop frequency response of /i controllers (top) and Skewed-^ 
controllers (bottom). 
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